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Abstract In this paper we present a grid of stellar evolution 
models, computed with an up-to-date physical description 
of the internal structure, using the Code d' Evolution Stel- 
laire Adaptatif et Modulaire (CESAM). The evolutionary 
sequences span from the pre-main sequence to the begin- 
ning of the Red Giant Branch and cover an interval of mass 
typical for low and intermediate mass stars. The chemical 
composition (both helium and metal abundance) is the solar 
one. The frequencies of oscillation, computed for specific 
stellar models of the grid using the Porto Oscillations Code 
(POSC), are also provided. 

This work was accomplished in order to support the 
preparation of the CoRoT mission within the Evolution and 
Seismic Tools Activity (CoRoT/ESTA). On the other hand, 
the grid can also be used, more generally, to interpret the 
observational properties of either individual stars or stel- 
lar populations. The grids (data and documentation) can be 
found at |http : / / www .astro . up . pt / cor ot /models 
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1 Introduction 

The computation of large grids of stellar evolution tracks 
is particularly useful to study a large number of stars for 
which photometric and/or spectroscopic observations are 
available. During the last 20 years, both the improve- 
ment of the physical description of stellar interiors (in par- 
ticularly, the opacities and the equation of state) and the 
computation facilities allowed many groups to produce a 
huge amou nt of stellar models. These cover from very-low 
mass stars Jchabrier and Baraffd lioooh to higher mass ones 



( Maeder and Mevnet 



2000), for different values of metal 



abundances and age s, including the solar standard models 
(IBahcall et al.l l2006). It is difficult to im agine that the lates t 
progress in stellar population synthesis (IGallart et alJl2005b 
or on the chemical evolution in the Galaxy would have bee n 
possible without the help of those grids (Pra ntzosl 120070 . 
^ c T > !fe a 5^ lication of a large number of stellar evolution mod- 
els to resolve single stars can, also, be very important in 
order to constrain the models themselves. For example, the 
comparison between isochrones and Hipparcos Population 
II stars, in the plan of the HR diagram, pointed out the 
limitation of tho se models to reprodu ce the observations of 
metal poor stars ( Lebreton et alJl2000h . 

With the advent of helio- and astero-seismology new ob- 
servational constraints for theoretical stellar models have 
become available. It is very well know how sensitive, for 
instance, the large and the small frequency separations are 
to both the global characteristics of stars (mass, radius and 
chemical composition) and to t he theoretical treatme nt of 
the stellar internal structure (e.g. lMonteiro et alJl2002l) . On 



the other hand, for fixed input physics, the models can be 
particularly interesting to predict the frequencies of oscil- 
lation of stars not (yet!) observed. The prediction can be 
important to select the stellar targets that will be measured. 
If the observational program is space based, those predic- 
tions are not only i mportant but absolu tely crucial to define 
the final target list ( Michel et al. 2006h . The grids reported 
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Table 1 Summa ry of the input physic s adopted to calculate the evolution models in the grids described in this work. Further details are 
also provided bv lLebreton et al. (2008). 



Item of the physics 


Selection 


Reference(s) 


Equation of State 


OPAL 


Rogers et al. (1996. 2005 Tables) 


Opacities 


OPAL 


Iglesias and Rogers (1996) + Alexander and Ferguson (1994) 


Nuclear reaction rates 


NACRE 


Anguloetal. ( 19991 


Convection 


MLT 


Bohm-Vitense (1958) + Henvev et al. (1965) 


Overshoot 


none 




Diffusion/settling 


none 




Mixture 


solar 


Grevesse and Noels (1993) 


Atmosphere 


gray 





here have already been us ed to support the seismic study of 



pre-m ain sequence stars dRipepi et alj l2007t iRuoppo et al 
2007h . 



In this work we briefly describe the grids of stellar 
models and their frequencies of oscillations produced in 
order to support the preparation of the CoRoT Mission 
(IFridlund et al.l 120061) . We start by identifying the physics 
and parameters that have been used to produce the models 
and then we discuss the data made available online. 



2 Parameters of the grids 

The present grid has been computed usi ng the cesam code 
(|Morellll997t iMorel and Lebretonll2008h . version 2k. In or- 
der to produce a grid with different initial conditions (birth- 
lin e) an extension o f the above version of CESAM, developed 
by iMarquesI (120081) . has been used. 



2.1 The input physics 

The input physics used to calculate the models is similar to 
the reference physic s adopted by the CoRoT/ESTA group 
dLebreton et al.ll2.008b . It is summarised in Table Q] The 
models consider neither microscopic diffusion of chemical 
elements nor overshoot of the convective core. 

2.2 Stellar parameters 

We have calculated evolution models with masses ranging 
from 0.8 to 8 Mq, from the contraction phase to the begin- 
ning of the red giant branch (RGB). All models were com- 
puted using the solar abundance of metals and helium and 
also the solar mixing length parameter (see Table [2] for a 
summary). The zero age main sequence (ZAMS) and the 
terminal age main sequence (TAMS) in the grids are defined 
as follows: 

ZAMS: when 99% of all energy is produce by nuclear re- 
actions; 



TAMS: when the abundance of hydrogen at the centre is 
X c = 0.01 ±0.0001. 



Table 2 Summary of the stellar parameters used to calculate the 
evolution models in the grids described in this work. Here Mq is 
the mass of the sun. 

Stellar parameter Value 

Helium abundance Y = 0.28 

Abundance of metals Z/X = 0.02857 

MLT parameter a = 1 .8 

Stellar mass 0.8 < M/Mq < 8 .0 



3 Initial conditions for stellar evolution 

Two options have been considered for the initial conditions 
of the evolution for each stellar mass. 

3.1 GRID A: evolution from fully convective spheres 

The initial models used for the pre-main sequence (PMS) 
evolution have tradicionaly been fully con vective models 
with arbitrarily large radii (about 50 Rq, e.g. lHavashil[l966l: 
Iben 1965). The reason behind this is that it was thought 



that stars were formed by the homologous contraction of the 
protostellar cloud; this non-hydrostatic phase was assumed 
to be so fast that the energy radiated was much lower than 
the gravitation al energy of th e cloud. A simple energy bud- 
get argument (IHavashil ll966) lead then to this initial large 
radius. 

This simple picture of star formation is n o longer held to 
be adequate dPalla and Stahlerlll99ll Il992h . Nevertheless, 
the high convenience of the so-called "classical" initial mod- 
els means that they continue to be widely used; besides, the 
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Fig. 1 HR diagrams for models of different masses at different 
ages: top panel - pre-main sequence evolution; bottom panel - 
main sequence evolution (full lines) and post-main sequence evo- 
lution (dotted lines). Only PMS models with an age above 1 Myr 
are shown. The number indicates the mass of the model in solar 
masses (Mq). 



influence of the initial conditions fades as the star evolves 
on the PMS. 

The example of the grid that uses this initial condition 
is illustrated in Fig. [T] in particular in the upper panel where 
the pre-main sequence evolution from an age of about 1 Myr 
is shown for models of different stellar masses. 

3.2 GRID B: evolution from a birthline 

It is now known that stars are formed in a very non- 
homologous way: an hydrostatic core is formed first, on 
to which matter from the parental cloud accretes, either di- 
rectly or through an accretion disk. Consequently, more re- 
alistic initial conditions are used for the calculation of this 
second grid. The most important reason to do so is that the 
PMS phase becomes shorter with increasing stellar mass, 
i.e. stars with higher masses are born nearer the ZAMS 



(see Fig . [2j ■ More realistic initial conditions are not relevant 
for stars with masses lower than ~ 1.5M©, but they become 
more important for stars with higher masses. In particular, 
stars with masses higher than ~ 1M® are born so near the 
ZAMS that they do not have a proper PMS phase. 



s 



TAMS 




Fig. 2 Stellar age at the ZAMS and TAMS for different stellar 
masses, in the two grids. The age at the ZAMS decreases must 
faster with stellar mass for the grid using the birthline (B) than 
the grid starting from the classical (politrope) initial condition (A). 
There is no significant difference for the age at the TAMS in the 
two grids. 

We started the calculation of the initial models with a 
0.5M© model evolved along the Hayashi track until its ra- 
dius was the same as t he radius of an accretin g protostar 
with the same mass in IPalla and Stahlerl (1199 ll) . Then ac- 
cretion is turned on at a constant rate of M = 10~ 5 M© yr _1 . 
The actual mass accretion rates vary with each individual 
case, and are not in any way constant JSchmeja and Klessenl 
2004). In particular, the mass accretion rate seams to in- 
crease with increasing mass. Given the goal of this grid, it 
is not necessary to include all the complicated phenomena 
of stellar birth in the calculation of the initial models; our 
accretion rate is intended to represent a typical value. 

The evolutionary track of the accreting protostar on the 
HR diagram is called the birthline. An initial model for PMS 
evolution is obtained when the accreting protostar reaches 
the mass of the model. Both, the birthline and some pre- 
main sequence evolutionary tracks, are illustrated in Fig. [3] 

An important influence on the location of the birthline 
is the deuterium content of the parental cloud. The mass- 
radius relation of the protostar is set to a high degree by 
deuterium burning. In "classical" PMS stars, the structural 
influence of deuterium burning is small: due to its low abun- 
dance, the deuterium burning phase is short. In accreting 
protostars, however, continuous accretion ensures that fresh 
deute rium is contin uously sup l ied (s ee IPalla and Stahlei 
199 lb . Following iLinskv et al.1 (|2006), we select a lower 
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4.0 
Log T„„ 



Fig. 3 HR diagram displaying the birthline (full line) and the pre- 
main sequence evolution of some models (the number indicates the 
mass of the model in solar masses M©). 



[D/H] number ratio than typically used; this was [D/H] = 
1 .5 x 1(T 5 , rather than 2.5 x 1(T 5 or 3 x 1(T 5 . 



4 Data 

The data that has been produced, as part of the grid, in- 
cludes: 



• evolutionary sequences (files of type * . datl and l* . HRl l. 

• stellar models at several ages (files of tvpe l* . OSCD . 



• oscillation frequencies (files of type * . f req). 
All files are available online for download at 

|http : / / www . astro . up . pt / cor ot /models/ cesam/| 

Specific information or data about the grids can also be 
provided on request from the authors. 

4. 1 Evolutionary tracks and stellar models 

The evolution of models of different stellar masses are pro- 
vided using the * . datl files (in same cases de l* . HRl output 
files from CESAM are also included). Here the evolution of 
the global parameters of the star is registered. The panels 
in Fig. 2] illustrate, as an example, how the size of the con- 
vective zones changes with time for models of 1 .0 M Q and 
3.0 Mq stars. The HR diagrams in Figs [T]|3] are also exam- 
ples of the data available in these files. 

For each value of stellar mass, several interior models 
at different ages can be found in the I* . OSCl files. For the 
details on the format and contents of the files please see 
the documentation available at the website. The files with 
the data on the internal structure of a star (with a specific 
mass) at a given time are grouped according to the follow- 
ing phases of evolution: 




Fig. 4 Evolution of the size of the convective regions (shaded 
area) with time, in models of GRID B, having masses M=1.0 Mm 
(top panel) and M=3M© (bottom panel). The full evolution from 
the birthline to the RGB are represented (dotted vertical lines rep- 
resent the ZAMS and the TAMS). 

prems - pre-main sequence evolution, 
ms - main sequence evolution, 
postms - post-main sequence evolution. 

Each model has a l* . OS C file containing all major quanti- 
ties describing the internal structure of the star (and used to 
calculate the seismic properties). 



4.2 Frequencies of oscillation 

For some of the models of the grid the freq uencies of oscil - 
lation were computed using the POSC code ( Monteirol2008 ). 
These are fhe |* . f req| files containing the basic mode pa- 
rameters. Only low degree modes are included (0 < I < 3) 
with mode orders from n — 1 and with a frequency below 
the acoustic cutoff frequency of the models. In same cases 
g-modes are also given (identified with a negative mode or- 
der). 

The boundary condition at the atmosphere used in the 
computations assumes an isot hermal atmosph ere at the top. 
For further details please see Mo nteirol ( 2008 ) and the doc- 
umentation available at the website. 



5 Concluding remarks 

In this paper we have reviewed the basic aspects of ref- 
erence grids of stellar models produced with the CESAM 
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code jMorel and Lebretoi3l2008h . Two grids are provided 
with different initial conditions: A - fully convective sphere 



(poli t rope); B - the birthline from Pall a and Stahleil dl991 



1992) . The remainin g aspects of the physics are fixed as in 



Lebreton et al.l (120081) . The grid also provides frequencies 
of oscillation for p- and g-modes of some of the models. 

The major goal of the grids is to support the prepara- 
tion and exploitation of space missions for asteroseismol- 
ogy and in particul ar the French CoRoT space mission 



jFridlund et al.ll2006h launched in December 2006. All data 
is available onlineH. 
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